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^ ■ Abstract 

O-f The two-dimensional free-boundary problem of steady periodic waves with vorticity is 

considered for water of finite depth. We investigate how a flow with small-amplitude Stokes 
waves on the free surface bifurcates from a horizontal parallel shear flow in which counter- 
currents may be present. The bifurcation mechanism is described in terms of a dispersion 
equation; namely, wavelengths of Stokes waves bifurcate from the values defined by the roots 
of this equation. The latter generalizes that for irrotational waves and involves only quantities 
given on the horizontal free surface of the initial parallel shear flow. Sufficient conditions 
guaranteeing the existence of roots of the dispersion equation are obtained. Two particular 
vorticity distributions are considered in order to illustrate general results. 

Mathematics Subject Classification (2010) 76B15 ■ 35Q35 

1 Introduction 

We study the two-dimensional nonlinear problem of steady waves in a horizontal open channel 
of uniform rectangular cross-section occupied by an inviscid, incompressible and heavy fluid, say, 
water. The water motion is assumed to be rotational which, according to observations, is the type of 
motion commonly occurring in nature (see, for example, [S51 [SI] and references cited therein) . Our 
aim is to consider the bifurcation mechanism resulting in the formation of Stokes waves (periodic 
waves whose profiles rise and fall exactly once per wavelength) on a horizontal free surface of a 
parallel shear flow in which counter-currents may be present. A detailed study of such flows is 
given in our previous paper 50^ under the assumptions that the flow velocity depends only on the 
vertical coordinate and it forms a basis for the present approach. Its crucial point is to determine 
a bifurcation wavelength for which purpose a dispersion equation is used. This technique is similar 
to the method that was earlier applied in our investigation of the irrotational Stokes waves (see 
[49] , section 8.2.1). To the authors' best knowledge, the dispersion equation derived here (see ((T8)) 
in section 1.4) was not used previously. In order to illustrate how the general results work for 
particular vorticity distributions several examples are presented. 



1.1 Statement of the problem 

Let a horizontal open channel of uniform rectangular cross-section be bounded below by the hor- 
izontal rigid bottom. Let water occupying the channel be bounded above by a free surface. The 
surface tension is neglected and the pressure is assumed to be constant on the free surface. The 
assumption that the water motion is two-dimensional and rotational and the incompressibility of 
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water allows us to seek the velocity field in the form (^'y,— ^'x), where {X,Y) are appropriate 
Cartesian coordinates and ^ is a stream function (see, for example, [52 and [55 for its definition). 
The vorticity distribution aj(^) is supposed to be a prescribed Cj^^^'^-function on IR, a S (0, 1), with 
bounded derivative. In the present paper, all variables are non-dimensional and chosen so that the 
constant volume rate of flow per unit span and the constant acceleration due to gravity are scaled 
to unity. For this purpose lengths and velocities are scaled to (Q'^/gY^'^ and {QgY^^^ respectively; 
here Q and g are the dimensional quantities for the rate of flow and the gravity acceleration, re- 
spectively. We recall that {Q^ / gY^^ is the depth of the critical uniform stream in the irrotational 
case (see, for example, [4]). 

Let {X, Y) be such that the X-axis lies in the longitudinal section of the canal's bottom and 
gravity acts in the negative F-direction. Moreover, let the frame of reference be chosen so that the 
velocity field and the unknown free surface are time-independent in these coordinates. Assuming 
that Y = S,{X) represents the free-surface profile (here ^ is a positive, continuous function), we 
denote by 2) the longitudinal section of the water domain, that is, 

D = {~oo <X < +00, 0<Y < CiX)}. 

Since the surface tension is neglected, ^ and ^ must satisfy the following free-boundary problem: 

*jfx + *yy+w(*) -0, iX,Y)eV; (1) 

«'(X,0) = 0, XeM; (2) 

^{X,aX)) = l, XeH; (3) 

\Vx,Y^{X,aX))\^ + 2aX) = 3r, X eU. (4) 

In the last relation (Bernoulli's equation), r > is the problem's parameter referred to as the 
non-dimensional total head (also known as Bernoulli's constant). This formulation has long been 
known (see, for example, |32j): its derivation from the governing equations and the assumptions 
about the boundary behaviour of water particles can be found in |16j . 

Considering the question of formation of Stokes waves on a flat free surface of a shear flow, we 
suppose that * e C2'"(I)) and ^ G C^-", < a < 1. 

1.2 Background 

Nonlinear theory of gravity water waves with vorticity has a long history which dates back to the 
1800s, when Gerstner [33] found his remarkable explicit solution (see [6] for a modern approach 
to this solution; in [36] some its properties are obtained). Notice that Scott Russell and Stokes 
published their pioneering works about irrotational water waves only in the 1840s (see [54] and [56], 
respectively). Nevertheless, rotational waves have been studied to a less extent than irrotational 
ones (see the recent survey paper [57], where an extensive bibliography is provided). Only a few 
articles treating waves with vorticity rigorously were published during the 20th century, and at 
least three of them are of lasting interest. As early as 1934, Dubreil-Jacotin |32j had proved the 
first existence result for these waves (in the 1950s, her work was extended by Goyon [34]). For 
this purpose she introduced a partial hodograph transform which is a very convenient tool to 
investigat rotational waves provided the horizontal component of the relative velocity in the flow 
does not change sign. In 1978, Keady and Norbury [45 obtained bounds on the total head of flow 
and free-surface profiles, but their results were proved only for a rather small class of vorticity 
distributions. In the authors' work [51', these bounds were generalized to the case of arbitrary 
Lipschitz distributions. 

During the past decade, a substantial body of rigorous results about waves with vorticity has 
appeared, the first of which was the article [16) by Constantin and Strauss. They used the partial 
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hodograph transform for obtaining a global branch of large-amplitude, Stokes waves. No counter- 
currents are allowed in this work as well as in the most of other papers that are briefly characterized 
by placing them into the following overlapping groups each of which covers a particular topic. 

• The existence of waves is shown through the local/global bifurcation mechanism in O [191 l20l 

[HSllSIllMlllIllMllMllSSllSniinil. 

• Regularity (in particular, analyticity) of the stream function and/or streamlines is investigated 

in [1312011311 Eg. 

• Symmetry of periodic and solitary waves is studied in [TUl [TTl [351 [Ml SO] ■ 

• Stability and instability of waves is considered in [TH] and [13] , respectively. 

• The works [Zl[ia[ni[Il[23[Ml[2a[2ll[30l[42[6Tl[Mll63[^ deal with some other properties 
of waves. 

• Deep-water waves are treated in [Tl] [25l [37l [41] . 

• Solitary waves are investigated in [35l [39l [40] . 

• The papers [H [21 [TH [2H1 [SS] are devoted to various forms of a dynamical system approach to 
waves with vorticity. 

• Numerical results about rotational waves are presented, for example, in [46} l47]. 

Now we turn to a few articles [20l [281 ISSl [30l [68], in which rotational waves on flows with 
counter-currents are studied (the last of these papers contains an extensive list of references and a 
brief review of the literature on water waves with vorticity). The works |20 [ I30 [ [68 ] are concerned 
with the case of a constant vorticity, whereas a linear vorticity is considered in [28l [29] . Since the 
results obtained in the cited papers are closely related to ours, we discuss them in greater detail. 

It should be emphasized that counter-currents are not considered as the principal point of those 
articles. Instead, the authors concentrate their attention on the appearance of internal stagnation 
points and the so-called critical layers. A critical layer is a connected subset of the water domain 
consisting of closed streamlines and stagnation points. See, for example, figures 1 and 3 in |68j 
and figure 3 in f20,. However, there exist degenerate cases like that shown in figure 2, [68 . From 
a geometric point of view critical layer is a horizontal array of cat's eye patterns anticipated by 
Kelvin [6^ (see also [22], p. 141). Therefore, such a layer, generally speaking, separates two other 
layers with opposite directions of flow. 

First, the existence of steady waves on a flow of constant vorticity in which a critical layer 
(and so, a counter-current) is present was established by Ehrnstrom and Villari [30], who studied 
streamlines and particle paths in the framework of linear theory. It was Wahlen [68' , who proved 
the existence of small-amplitude Stokes waves with a constant vorticity and discovered that in the 
reference frame moving with the wave there is a critical layer. (As is discussed in [68], the term 
itself has a long history, but was mainly used in the framework of mathematical models other 
than rotational water waves.) Recently, Constantin and Varvaruca [20) proposed an approach to 
Stokes waves with a constant vorticity that differs from that in |68] and is based on the conformal 
mapping technique that imposes no restriction on the geometry of the free surface profile which, 
in particular, can be overhanging. This technique, as the authors claim, opens up the way to using 
global bifurcation theory. Besides, the physical relevance of the problem with a constant vorticity 
is discussed in [20] . 

In the articles [28l [29], Ehrnstrom et al. studied the case of linear vorticity when multiple 
counter-currents and critical layers exist. In [29], the results obtained in |68j were extended, yielding 
the existence of small-amplitude waves with arbitrarily many critical layers (including the so-called 
bichromatic waves). Qualitative and some quantitative properties of waves the existence of which 
is proved in [29] are investigated in [28]. In particular, a classification of vorticity distributions is 
proposed and the bifurcation relations are given for different cases. 
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1.3 Stream solutions and their basic properties 

A pair {u{Y), h), where h = const, is called a stream solution of problem (Il|)-(j4]) when 

= u{Y) and ^{X) = h 

satisfy this problem. Such solutions are studied in detail in [50], where they are also discussed in 
terms of the unified theory of conjugate fiows developed by Benjamin [3]; here we restrict ourselves 
to a summary of those results. 

Seeking stream solutions, we write problem (UJ-lHl) in the form: 

u" + w(u) = 0, ii(0) = 0, u{h)^l, \u'{h)\^ + 2h^3r, (5) 

where u' — uy- By U{Y;s), we denote a unique solution existing on the whole IR for the first 
equation ([5]) complemented by the following Cauchy data: 

C/(0; s) = 0, C/'(0; s) = s, where s e R. (6) 

It is reasonable to use s for parametrising solutions, and it occurs that all stream solutions are 
parameterised by 

s> So = ^2 ma.x^n{T) , where 17 (r) = J uj{t)dt. (7) 

Integrating the first equation ([5]) for U{Y; s), one obtains 

([/')' + 2 r!(C/) = s2. (8) 

Moreover, we get that if h{s) > is found for some s so that 

U{h{s);s) = l and s^ - 2n{U{h{s); s)) +2h{s) ^ 3r (9) 

hold simultaneously, then the pair {U{Y;s), h{s)) is a stream solution (note that in the second 
relation the expression from ^ is substituted). These solutions exist only for r > rc, where 
Tc is defined by formula (|15p below, and the set of stream solutions corresponding to a particular 
r can be obtained by virtue of the following procedure (see details in [50]). 

One takes some s > sq [see ([7])] and finds U{Y; s) from the Cauchy problem described above. 
In fact, s depends on r and we shall describe below how to get s for a given r (note that several 
values of s might correspond to some particular r). Then one obtains the values of depth for 
which purpose the first relation ^ serves. Let us outline the corresponding scheme involving two 
auxiliary quantities which we denote t±{s) and y±{s) (see |50| . section 3, for their properties). 

If the equation 2r2(T) — with s > sq has finite positive and negative roots, then by t+(s) 
[t_(s)] we denote the smallest positive root [the largest negative root, respectively]. If there is no 
finite positive [negative] root, then we put T4-(s) — +oo [t_(s) = — oo, respectively]. Notice that 
''"+(s) > 1 if and only if s > sq. Now we set 

V±(s) — / — — for S > Sq. 

^ Jo ^s^ - 2Q{t) 

The introduced quantities are such that (?/_(s), y+(s)) is the maximal interval, where U{Y; s) in- 
creases strictly monotonically, and t+(s) [t_ (s)] is the supremum [infimum, respectively] of U (Y; s) 
on this interval. More precisely, U (Y; s) has the following properties. 

• If y+{s) = +00 and y-{s) = — oo, then U{Y;s) increases strictly monotonically for all Y E H. 
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• If 2/-(s) = — oo and y+{s) < +00, then U{Y; s) is bimonotonic and attains its maximum t_|_(s) 
at F = y+is). 

• If y-{s) > —00 and y+(s) = +00, then U{Y;s) is bimonotonic and attains its minimum t_(s) 
atr = y_(s). 

• If both y+{s) and y-{s) are finite, then s) is harmonic-hke; it attains one of its minima at 
Y = y-{s) and one of its maxima atY = y+{s). Moreover, U{Y; s) increases strictly monotonicahy 
from r_(s) to t+(s) on [y_(s), y+{s)]. 

Furthermore, it is proved in [50) . section 4.2, that for s > sq 

dr 

is the smahest positive value of the depth such that U{h{s); s) = 1. By /i^^''(s), j = 0, 1, . . . , we 
denote the whole sequence (possibly finite) of depth values such that U (^/ij^"* (s); = 1 for a given 
s, and we have 

h^^,\s)^h{s) + 2k[y+is)-y^{s)], j ^ 2k; (11) 
hiV+i{s)^h{s) + 2[y+{s)-h{s)]+2k[y+{s)-y^{s)], j - 2fc + 1. (12) 

Here, fc = 0, 1, . . . , and so h{s) [see formula ([TU|] is included into (|lll) . where it appears as h'^\s). 
If both y+{s) and j/_ (s) are finite, then formulae pj|) and give finite values for all fc = 0, 1, . . . . 
Otherwise the first of them gives a finite value only for fc = 0. Formula also gives a finite 
value for = provided y+(s) is finite. Note that both formulae coincide when y+{s) = h{s) that 
is equivalent to the equality t+(,s) = 1. 

For s > So, there exists the following sequence of solutions 

h\~\s)^hl+\s)~2y^{s), J =0,1,..., (13) 

of the equation U{h{s); — s) = 1, and so this sequence also gives flow depths. The number of finite 
elements in the sequence (jl3p depends on whether j/+(s) and y-{s) arc finite or not. It is clear 
that the value /i^ \s) is infinite when either h'j^\s) ~ +00 or y-{s) = —00, whereas /ij \s) is 
finite otherwise. 

The subscript (superscript) at h^j^\s) points at the number of counter-currents (the direction 
of fiow in the near-bottom layer, respectively) in the following way. If the depth is equal to 
then the number of layers with alternating directions of flow is equal to j -|- 1, whereas the direction 
of fiow in the near-bottom layer coincides with the direction of the so-called critical unidirectional 
flow, which always exists, corresponds to r = Tc and has the minimum depth among all possible 
flows. If the depth is equal to /i^ \ then there are j + 2 layers with alternating directions of flow 
and the near-bottom layer has the direction opposite to that of the critical flow. 

Let us summarize our algorithm for finding stream solutions. We find appropriate values of s 
from the second relation ^ (Bernoulli's equation): 

3^(±)(s)=r, where Ji'f^s) = \\s^- 2 + 2 hf\s)], j=0,l,..., (14) 
and the last expression is obtained with by virtue of ([5]). Then the stream solution 

{u{Y;s),h^+\s)) [{u{Y;~s),h^-\s)y_ 
corresponds to the root s of equation (fT4)) with the subscript j and the superscript -I- [- 
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Let us turn to the definition of rc- Since for every s > sq the sequences 3l^^^\s) and \s) are 
increasing (they might be infinite or finite depending on the vorticity distribution), we have that 
the graph of 3?[,+^ (s) (it is convex as is shown in |50) . section 5.1) lies below those corresponding 
to other functions 3?j''''(s). Therefore, we put 

rc = min3l{i+^(s) (15) 

s>so 

and this value is always attained at some Sc > sq. Moreover, for all r > rc some equations in the 
sequence ((T4)l have roots; at least two such roots exist for every r. It was said above that for each 
s obtained from the sequence of equations (ITil) one finds U{Y;s) and h{s). Thus, the set of all 
stream solutions corresponding to r > rc is determined. 
Finally, it must be mentioned that if 

ro= lim 3l^+^(s) (16) 

is finite, then among the flows corresponding to every r > tq there exists at least one flow that 
has a counter-current. These flows are defined by stream solutions obtained above and flows with 
several counter-currents can exist among them. 



1.4 Main results 

The existence theorem for Stokes waves involves three assumptions. Our first assumption concerns 
determining a shear flow from which Stokes waves bifurcate, namely, we suppose the following. 

(I) The inequality r > rc holds. As was said above, it guarantees that some equations in the 
sequence ([T4|) have roots. Moreover, at least two roots exist for every r > rc, and we choose those 
which are greater than sq. Let s* be one of them, and so we have that U'{h; s*) ^ (see [50]). 

The next assumption is related to the so-called dispersion equation and properties of its roots. 
Similarly to the case of irrotational waves every of these roots determines a bifurcation wavelength 
as a function of r and h (of course, h itself depends on r). In the case of waves with vorticity, 
more than one bifurcation wavelength can exist for a given choice of a stream solution and this 
depends on the vorticity distribution (see section 5). On the contrary, only one wavelength exists 
for a given admissible r in the irrotational case because the corresponding dispersion equation can 
be written as follows (see [49], p. 478): 

r coth hr = . (17) 

Here /i > 1 is defined by r and its value is equal to the non-dimensional depth of the subcritical 
uniform stream (for convenience, the depth of the critical stream is scaled to unity). It is clear 
that (|17p has only one positive root which is equal to the wavenumber 2 7r/Ao corresponding to 
the bifurcation wavelength. It is worth mentioning that the problem about periodic irrotational 
waves has another statement that leads to a different form of the dispersion equation. The latter 
involves the phase velocity of waves (see, for example, formula (1) in [69] . ch. 71); however, this 
form is less convenient for generalising to the case of waves with vorticity. 

In the rotational case, the stream solution {U (•; s*), /i(s*)) corresponding to the chosen yields 
the following dispersion equation (see the derivation in section 2.4): 

cr(r)=0, where a{T) ^ K-f'{h,T) - (18) 

In the last function, k = U'(h), whereas 7(y,T) solves the following problem: 

-7"-|-[t2-c^'([/)]7 = 0, 7(0,r) = 0, jih,T)^l, (19) 
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in which h depends on r through the root of Bernoulli's equation. Note also that 



K = ±VZ^^^:^ when ± U' {h) > 0. 

If is not a Dirichlet eigenvalue of the operator d^/dY^ + w'(?7) (there might exist a finite 
number of such eigenvalues), then problem ()19p has a unique solution. Therefore, 7(y, r) and (T(r) 
are defined for all these values of t and are smooth and even functions of r. There are vorticity 
distributions for which the Dirichlet spectrum is not empty, and so ^(Y, r) and (t(t) are not defined 
for some values of r. For example, this is true when w(t) is a linear function positive for r > 
(see section 5.2). 

Let us list some properties of equation (fT8|) . First, the function (t(t) is well-defined only when 
s* 7^ So because otherwise k vanishes. (Note that the equality = sq, excluded by assumption 

(I) , holds only when tq defined by formula (|16p is finite.) Second, if lj vanishes identically, then 
(fTSj) coincides with (|17p because k = h^^ and U ~ Y/h in this case. Then the solution of (|19p is 
7 = sinhFr/ sinh /ir which gives (fT7|). Third, if zero is not a Dirichlet eigenvalue of d"^ / dY^ +uj' {U) 
(this is equivalent to the inequality (dh/ds) (s*) ^0 and follows from assumption (III) formulated 
below), then the following equality holds (see formula ([55|) in Proposition 3.3): 

2 \d5l fdh - 
<7(0) = -— — (s*) / -^-(s*) ■ 
3k [ ds / ds 

For the sake of brevity, here and below Jl stands for the function ^f^"^ that defines the stream 
solution {U,h) which is involved in (|T8| and p9)) . The last equality is used for investigating 
whether equation (|18p has roots. In particular, the existence of a positive root in the absence 
of the above mentioned Dirichlet eigenvalues follows from positivity of the fraction in the square 
brackets. On the other hand, if there is a Dirichlet eigenvalue, then p8|) has at least one root 
greater than this eigenvalue [see Proposition 3.3 (ii)]. Moreover, if the fraction in the square 
brackets is positive, then there also exists a root between zero and this eigenvalue. Finally, for 
any r > rc the equation 3?g^''(s) = r has a root > Sc- If the stream solution ([/, h) corresponds 
to this s* and a defined by ([/, h) is a continuous function, then equation ([TS]) has no positive 
solutions. The last fact has the well-known analogue for zero vorticity: only solitary waves exist 
in the supercritical case (see, for example, f2| and J48l). 

Now we turn to the second assumption. 

(II) The dispersion equation (|18p has at least one positive root, say, tq such that: (a) (t'{tq) ^ 0, 
which means that this root is simple; (b) none of the values r^jk (fc = 1, 2, . . . ) is a root of (fTS]) . 

In section 3, we consider conditions guaranteeing that this assumption is fulfilled. Now we are 
in a position to formulate the following. 

Definition of a bifurcation wavelength. Let assumptions (I) and (II) hold, and let tq be a 
root of (fT8| existing by assumption (II). Then we put Aq = 27r/To, and so Aq is a wavelength of 
linear rotational waves on the shear flow described by the stream solution ([/, K) corresponding to 
s* introduced after assumption (I). 

If equation (fT5|) has more than one root satisfying conditions (a) and (b) of assumption (II), 
then each of these roots defines a wavelength of linear waves with vorticity and all these waves 
exist on the free surface of one and the same shear flow. Our last assumption is as follows. 

(III) If the boundary value problem 

V5xx + </3yy +t^'(t/)'^ = inRx (0,/i), v'(X,0)==0, ip{X,h') = Q, 
has a weak solution such that it is Ao-periodic in A", then this solution is trivial. 
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In particular, this implies that zero is not a Dirichlet eigenvalue of the operator d^/dy^ + io'{U) 
considered on {0,h), and so the function cr(T) is defined for t belonging to a neighbourhood of 
zero. 

Let us describe function spaces that are used in what follows. For a non-negative integer k and 
a G (0, 1) by H^'" we denote the space of Ag-periodic, even C*^'" -functions on R (cf. 49 , p. 478). 
Furthermore, let the strip D be a C'^'^-domain whose boundary is Ao-periodic and even in X. 
Then n^'"(2)) is the space of C'^'^-functions on D which are A-periodic and even in X. Now we 
are in a position to formulate the following result. 

Main Theorem. Let oj G C^'"(]R), a e (0, 1), and let r > satisfy assumption (I) and be such 
that assumptions (II)~(III) hold for a stream solution {U, h) defined by r. Then there exist e > 
and a continuous mapping 1 1— >■ (A(i), C('i 0) ffom {t £'Si:\t\ £ (0, e)} to a neighbourhood of (0, 0) 



II X n^'", where Aq = I'k jr^ and tq is a positive root of equation (1181) existing by assumption 



(II) . Moreover, the components of this mapping have the following properties. 

(i) The first component X{t) ^ and X(t) as \t\ — > 0, whereas Ci^^'t) ^ for every t and 
X G [— Ao/2, Ao/2], and the asymptotic formula holds 

'2ilTX 

(:{x,t)=tcos— + C*ix,t), (20) 

where ||C*(-, t)||j-[2,a = o{t) as \t\ — > 0. 

(ii) // the upper boundary of D has the form 

'''' 

that is, t) 0, A-periodic function of X with A — Aq[1 + X(t)], then problem ((T])-([3|) has a 

solution '^{■,t) G n^'"('I)) such that the pair {^{X,Y,t), ^{X,t)) satisfies condition Moreover, 

^i,{X,Y,t) = u(Y^^^+^4X,Y,t), (22) 

where ||**(-, t)!!^^,^^^,) = 0{t) ast^O. 

This theorem means that (vf, 5) is a A-periodic Stokes- wave solution of problem ([I])-©; for- 
mulae (pn|) and (|22p give the asymptotic representations for its components. Thus, this solution 
is a perturbation of the stream solution ([/, h). Relations ([20)) and (I?!]) immediately yield this 
for ^, whereas for this fact will be shown in the theorem's proof, for which purpose the form of 
function ^I^* will be specified (see section 4.2). Moreover, ^) belongs to a family of solutions 
describing Stokes waves with wavelengths close to Aq provided the latter is defined by a simple 
root of the dispersion equation ^TE\\ . Note that if is finite and r > tq, then (U, h) describes a 
shear fiow with counter-currents (see (50t and section 5, where examples are considered). 

In conclusion of this section, it is worth emphasizing that the mechanism of wave bifurcation 
is substantially more complicated for waves with vorticity comparing with that for irrotational 
waves. First, for every r > rc there are exactly two uniform fiows (sub- and supercritical) in the 
irrotational case, whereas there can be as many rotational fiows of constant depth as one pleases 
which, for example, is the case for linear positive vorticity (see section 5.2 below). Second, for a 
given depth of the subcritical irrotational uniform fiow the dispersion equation ()17p defines only 
a single value of the bifurcation wavelength. On the other hand, the structure of roots of the 
rotational dispersion equation (jlSp depends on the value of depth (the second component of the 
stream solution), which, in its turn, depends on r > rc. This results in the fact that the number 
of values of the bifurcation wavelength can be as large as one pleases, which again is the case for 
linear positive vorticity. 
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2 Operator form of problem ([T])— (jlj); dispersion equation 



We prove Main Theorem using the Crandall-Rabinowitz theorem (see Theorem 4.1 in section 4). 
In order to apply it to problem ([IJ-Q, the latter must be transformed into a nonlinear operator 
equation. 

2.1 Reformulation of the problem 

First, we have to choose a horizontal shear flow from which Stokes waves can bifurcate. According 
to assumption (I), we fix r > Tc, where the critical value is defined by formula (1151) and depends 
on the vorticity distribution. Then equation ()14|) has a solution s* at least for one function ^^j^"*, 
J = 0, 1, ... , and the same assumption implies that 7^ so- The stream solution corresponding to 
s* is as follows: 

{U^^\Y-s.),hf\s.)). (23) 

Here the second component is given by one of the formulae ([Tl]) -(fT3 |) and expressions for the first 
component can be found in j50j , section 3, whereas its properties are described above in section 1.3. 
For the sake of brevity, this solution will be denoted (C/, h) in what follows, whereas 3? will stand 
for the function 3?^*-' defining {U,h). Note that U'{h,s.i) ^ 0, which is also a consequence of 

assumption (I). If for some r > equation has roots for several functions 3?j-*\ then each 
of these roots can be used for obtaining a dispersion equation through the corresponding solution 
(j23[) . If assumption (II) is fulfilled for roots of several dispersion equations, then more than one 
bifurcation wavelength corresponds to the chosen r. 

Now we reformulate the problem in order to introduce a small parameter so that the Crandall- 
Rabinowitz theorem is applicable to a single operator equation resulting from the amended state- 
ment. Notice that the chosen stream solution (?7, h) defines the wavelength Aq > through the 
dispersion equation (HH]), because {U,h) is taken so that assumption (II) holds for ([TS| and its 
roots. Let A denote the wavelength of the sought Stokes waves that perturb the free surface Y — h 
of the shear flow. Assuming that A is close to Aq, we put A = Ao(l + A), and so A is a small 
parameter to be found along with the stream function and the wave profile. 

For transforming problem (Il])-(l4]) into a new one depending on A, we introduce the following 
variables: 

- = y=^, ^(a.) = ^^, ^{x,y) = nX,Y). (24) 

Then D is changed to the curved strip D — {~oo < x < +00, < y < ri{x)}, and the problem 
takes the form: 

V'.. + V'yy + (l + A)2c^(V) = 0, {x,y)eD; (25) 
'i{j{x,0)^0, xeIR; (26) 
7/;(a;,r/(x)) = 1, xelR; (27) 
\Vr,yijix, y)p + 2(1 + X)^y = 3r(l + A)^, y = r;(x), xelR. (28) 

Now we seek (ipyr]) and A, so that rj is a non-constant, Ao-periodic, even function, whereas ip is 
Ao-periodic and even in x. 

In view of the boundary condition (P7)) . equation ([^5]) takes the form: 

$^(x,?7(x))-(l + A)[3r-2(l + A)ry(a;)]^/'=0, xeR, (29) 
on 

provided U'{h) > 0; here n is the unit normal directed outwards of D. If U'{h) < 0, then we have 
^(x,?7(x)) + (l + A)[3r-2(l + A)ry(:E)]i/2 = 0, xeR (30) 
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instead of (P^ . In what follows, we restrict our considerations to (I^Hl) . just formulating the results 
that concern pop . 



2.2 Reduction to a single operator equation 



As in the case of irrotational waves (see [49], section 3), it is convenient to reduce problem (|25l) - ([28)) 
to a single equation. The first step is to transform the curved strip D into S* = II x (0, ft.). For 
this purpose we change the vertical coordinate y to the following one: 



h 



y 



r]{x) 



or z — Y— — - according to (IMl) 



(31) 



Setting u{z) = U{Y), we define new unknown functions C{x) on R and (t>{x, z) on S as follows 

h 



C{x) = r]{x) 



1 + A 



and 



{x,z) ^u{z)-'ijj(^x,^'q{x)j 



(32) 



Thus C and 4> are Ao-periodic and even in x functions, but they are small when a perturbed flow 
is close to that defined by the pair (u, /i), which is nothing else than the stream solution {U,h) 
written in the variables (x, z). 

The next step is to describe how (j) depends on A and which, in its turn, defines the mapping 
(A,^) I— >■ (j)z{x,K). Therefore, the latter function is considered in the final equations which follow 
from ((29)) and pO|) as the image of an operator applied to C,. We begin with the weak setting of 
problem (P5)) - ([77|) . and the corresponding integral identity is as follows: 



[i'xVx + '4>yVy - (1 + A)^a;(V') v] dx dy = 0. 



(33) 



It is obtained from equation (P5|) by integration by parts and is valid for all v G W^^^{'D) vanishing 
on (91); by W^^^{'D) we denote the space each element of which belongs to W^'^{K) for some 
bounded open subset K C D. 

In identity psp . we change the variables (x, y) and the function tp to (x, z) and 0, respectively. 
Moreover, we take into account the Ao-periodicity and evenness of rj, and the fact that (j) is also 
periodic and even function of x. This allows us to integrate over (— Ao/2, Ao/2) instead of IR (it is 
convenient to take (—Ao/2, Ao/2) as an interval of periodicity), and so we arrive at the identity 



Ao/2 rh 
-Ao/2 Jo 



Dx [Uz + (Pz) 

V 



zrjx 



-{1 + \)'^LU{U + (p)v 



+ ( - ) iUz+ (j}z)Vz 



rj dx dz 



= 0. 



(34) 



It is valid for all v G VF^'^(5) that are vanishing on dS and Ao-periodic in x. Furthermore, we 
have that 



Ao/2 



Ao/2 Jo 



h dx dz 



UzVz 



Ao/2 
-Ao/2 Jo 



^ , . hdxdz 

UJ(U) V , 



because u in p4p is the chosen stream solution. Using this equality and the fact that rjx = Cxy we 
transform (134^ into the following form: 

X - —{Uz + (j^z) 

V 






zCx 






Vx Vz 











- uj{u) - (1 + X) u}(u + , 



rjdxdz 



= 0. 
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Let us rewrite this integral identity as follows: 



zC,x 



zC,x \ f h^, 

V ) \n, 



-(1 + Xf-\jj{u + 0) - uj{u)\v 



zCx 



(l + A)3[2/.+ (l + A)C]C , . \vdxdz 

5 —U}[U)V} ; . (35) 



Here v belongs to the same class as in (p4|. whereas 77 is the shorthand for the first term in the 
right-hand side in the first formula (|32|) (note that it depends on ( and A). 

The existence of cj) satisfying ([55]) is a consequence of assumption (III) provided C and A are 
sufficiently small. Precisely, this is formulated as the second part of the following assertion proved 
at the end of section 2.3. 

Proposition 2.1. Let assumption (III) hold, and let the pair (A, C) belongs to a sufficiently small 
neighbourhood of (0, 0) m IR x 11^'". If f ^ 11^'" (S*), then the weak problem 



zQx 



zCx 



(f>zVz - (1 + A)2 [uj{u + (f>)~ uj{u)] V 



•q dx Az 



f V Axdz . 



(36) 



where v and rj are the same as in psp . has the (S) -solution cj) that satisfies the estimate 
\\(j)\\jj2,o,,g, < C ||/||]-[0,Q (5), where the constant C is independent of f. 

Aq y > Ag 

In particular, the mapping (A, C) > defined by (I35p continuously maps a small neighbourhood 
of (0,0) inlRx n^'" to the subspace ofIl^^{S) whose elements vanish on dS; again the estimate 
\\4>\\„2,c, ,g. < C ||C||j-[2,c« holds with C independent ofC- 

Aq \ / Aq 

Now we turn to deriving a single operator equation for ( from the Bernoulli's equation 
Applying ([5^ . it is straightforward to get 

Furthermore, (j)x{x, h) vanishes identically, and so ([29]) takes the form: 

h 



[uz{h) + (j)z{x, h)] 



Tjix) 



(l + A)[3r-2(l + A)ry(a;)]^/^ = 0, xeU. 



As above 77 stands for the first term in the right-hand side of the first formula (|32p . Finally, we 
put K = Uz{h) [it is equal to (3r — 2/i)^/^ when Uz{h) > 0], thus obtaining after simple algebra the 
following form of Bernoulli's equation: 



(j)z {x,h) 



[k2-2(1 + A)C(^)]'^' + k = 



(37) 



Here x e (— Ao/2, Ao/2), the values of h and k are given [of course, both of them are functions of 
r because h depends on r through the root s* of equation (fT4|) ]. whereas C and A are unknown. 
Moreover, (f) depends on C and A through the weak formulation (1551) . and so ^^(a;, h) is the image 
of a nonlinear operator mapping (C, A) into this function on (—Ao/2. Ao/2). 
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We recall that equation (P7)) is obtained under the assumption that U'{h) > 0. Supposing that 
U'{h) < 0, we have n = Uz{h) ~ — (3r — 2ft,) and so Bernoulli's equation takes the following 
form 

[K2-2(l + A)C(a;)]'^' + K = 0, (38) 
instead of ([37|: here again x G (— Ao/2, Ao/2). 



4>zix, h) 















V 


da; dz = — 


r 1 









2.3 Linearization of equations ( I37p and ( 1381) 

Here we linearize equations (I37p and (I38p . thus arriving to a single linear equation whose role is 
twofold. First, it serves for deriving the dispersion equation ([T5)) (see section 2.4); second, it is used 
when one applies the Crandall-Rabinowitz theorem (see Theorem 4.1 in section 4) to equations 
([571) and ([55]). 

Proposition 2.2. T/ie linearized form of equations p7p anc? psp is as follows: 

fi^\x,h)- (j--^^C{x)^0, a; e (-Ao/2, Ao/2). (39) 

_ffere lya*-^-* — the leading term in the expansion of 4> — is an even, A^-periodic function of x that 
satisfies the integral identity: 

h 

[zUzCxVx + 2uj{u)(v]dxdz, (40) 

It j-fio Jo 

where v belongs to the same class as in (j34p . 

Proof. The fact that (ys'-^' is even and Ao-periodic in a; is a direct consequence of (|39p because ( 
has these properties. 

Now we turn to deriving ((5^ . In (P7)) and ([35]), the only linear term with respect to (j) is the 
derivative (j>z{x, h). Besides, it is straightforward to check that (l40l) is obtained by linearization of 
the integral identity psp defining (f). Thus, 4>z{x, h) must be changed to ip^z\x, h) in the linearized 
equation. 

Expanding the second term in ([37l) and ([38)) . we see that k cancels [one has to take into account 
that ~ —K in (j38p ]. whereas the first- and second-order terms containing A, C and (^^ are 

respectively. Now we change the second term in ([37]) and (|55)) to the first of these expressions; the 
first term in these equations we change to ipi^\x,h). Thus we arrive at the linearized equation 



Remark 2.2. Using assumption (III) and the well-known results on the regularity of solutions to 
elliptic boundary value problems (see, for example, [1]), one obtains that the problem 

V'xx + Vzz+uj'iu)ip ^ f in (p{x,0)=0, ip{x,h)=0, 

has a unique solution that belongs to 11^'" (S') provided / G 11°'" (S*). Moreover, the estimate 

ll^lln-(5)^^ll/lln-(5) (42) 



holds with the constant C independent of /. 
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Of course, this assertion is true for the following particular case: 

'^i'i+'^iV+^'M'^''^ = /i"'[^".C..-2w(u)C] in^, (^W(^,0) = 0, = 0, (43) 

where C is an even, Ao-periodic C^'"-function. However, the norm in the right-hand side of inequal- 
ity (1321) must be changed to ||C|ln^'° ^^^^ case. Note that the latter problem is the differential 

Aq 

form of (gOl)- 

Proof of Proposition 2.1. Let us consider the nonlinear and linear operators defined by the left- 
hand sides of the integral identities ([36]) and (|40)) . respectively. Since the pair {XX) belongs to a 
sufficiently small neighbourhood of (0, 0) in H x H^'", it is easy to see that the nonlinear operator 
is a small perturbation of the linear one. According to assumption (III), the linear operator 
is invertible (see Remark 2.2). Combining this fact and the contraction mapping theorem, one 
obtains the first assertion of proposition. Moreover, the estimate of |j(/)|jj-[2,c. .a-i is a consequence of 



inequality (j42|) . 

Furthermore, we use integration by parts in order to write the right-hand side of (|35p in the 
same form as the right-hand side of ([55)1 . Then the second assertion of proposition follows by 
applying the first assertion to the weak problem ([35]). 



2.4 Derivation of the dispersion equation (1181) 

In order to obtain the dispersion equation for waves with vorticity we apply the same procedure 
as in the case of irrotational waves (see [15], pp. 466 and 478, and [SS], PP- 448-449); that is, 
we use the linearized problem. In the present framework, the original problem is equivalently 
reduced either to equation ([37]) or to ([38| (which of these equations arises depends on the vorticity 
distribution). Both of them give the same linear relation (155]) between ( and ip^^^ (the latter is the 
leading term in the expansion of (f>). Along with (|39|). both these functions appear in the boundary 
value problem ([43]) . Hence for deriving equation ([T8]) it is sufficient to exclude them from ([39]) and 
(03]). 

Since C and ip^^^ are even functions of x, we seek them in the form similar to the principal 
term in formula ([20]) . that is, ({t)costx and (p^-^^r, z) cos tx, respectively. Then we get that the 
following relation holds for r G R: 

0{}''{T,h) — ( — j C{t): where (p''^''{T,z) satisfies the problem (44) 

-Pi\^ + [t^ -uj'{u)]ip^^^ ^h-^[T^zu, + 2uj{u)]C, ze{0,h); p^^^ {t,0) = 0, ^(^^ (r, /i) = 0. (45) 

According to assumption (HI), this problem has a unique solution. Let G{z,t',T) denote Green's 
function of problem (|45p . then we have that 



^(1) z) = ^ r G{z, t; t) [t2 tut + 2 oj{u{t))\ dt. 



Let us show that 

G.{z,t-r)\^^,^^-l{t,T), (46) 

where 7 is a solution of problem ([19]) . Denoting by L the differential operator in the left-hand side 
of equation in problem (|19p [the same operator appears in the equation of problem (|45p ] , we write 
the following representation: 

7(t,T)= / [7(z,T)(iG)(z,i;T)-G(z,t;T)(i7)(z,T)]dz. 
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It is clear that the integral is equal to 

[lz{.z,T)G{z,t;T) - Gz{z,t;T)-f{z,T)]l~^^ ^ -6*2(2, i; r)|^^^. 

The latter equality is a consequence of the boundary conditions (fT9)) and (|45|l . Combining the last 
two formulae, we arrive at 

Now we differentiate the formula expressing (p^^^r^z) with respect to z and use (j46p . thus 
obtaining 

(^ii)(T,/i) = -ao(T)C(r), where (To (r) = /i^W 7(2, r) [r^zu^ + 2 w(w(z))] dz. (47) 



Comparing this and ()44p . we get the following form of the dispersion equation: 

K 1 
f^o(T) + --- = 0. 

In order to obtain ([T5|. it remains to show that the left-hand side is equal to 

ct(t) = k7'(/i,t) - +w(l). (48) 

Problem and the problem for u allow us to calculate the integral that defines cro(T). Since 
uj{u) = —Uzz on the interval (0, ft.), we have that UzUj'{u) — —Uzzz there. Besides, the first relation 
gives that r^7 = 7^2 + ^ijj'{u). From these facts we get that 



o-o(t) = -y {z [7^2 + 7w'(u)] - 2 7M22}dz = - y {zl-yzzUz - ^ Uzzz] - '^1 Uzz} dz. 
Integrating by parts and applying the boundary conditions for 7, we obtain 

ZJUzzzdz ^ -huj{l) - {zj)zUzzdz. 

Jo 

Here we also used the relations Uzz{h) = —oj{u{h)) = — a;(l), which follow from the definition of 
u. This gives 

If'' If'' 

(To(t) = W(l) + - y [z{jzzUz + ■^zUzz) - JUzz]dz ^ Uj{l) + - J [z (^zUz) ^ - 1 Uzz\dz. 

Integrating by parts in the last integral, we see that it is equal to the integrated terms only; 
furthermore, taking into account the boundary conditions for 7 and the definition of k, we arrive 
at formula 



3 On roots of the dispersion equation (1181) 

In this section, we consider the dispersion equation (jl8l) corresponding to some stream solution. 
First, we prove sufficient conditions which guarantee that (|18p has at least one positive root, and 
so the first part of assumption (II) is fulfilled. Second, we show that ([T5)) has no roots in the case 
when this equation corresponds to a stream solution defined by a supercritical value of s, that is, 
s is greater than Sc and satisfies the equation 3?g~''-'(s) = r, where r is any number greater than 
r^. The latter fact is analogous to what is well known for zero vorticity, namely, that only solitary 
waves exist in the supercritical case (see, for example, j2] and [48]). 
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3.1 Sufficient conditions for solubility of equation ( llSp 

Let assumption (I) hold, and so there exists that solves equation (|14p for some 31, and let {U, h) 
be the stream solution corresponding to s* . (Here we use the simplified notation introduced above, 
namely, {U, h) stands for the stream solution, whereas 3? denotes the left-hand side of equation 
(fT4|) whose root s* defines this solution.) We begin with two auxiliary assertions concerning the 
function a defined by ([/, h). The first of them describes the behaviour of cr(T) at infinity, and the 
second assertion gives the asymptotics of j{Y, t) near an eigenvalue (under the assumption that 
it exists) of the operator d^/dY^ +uj'{U) with the Dirichlet boundary conditions. (We recall that 
only a finite number of such eigenvalues can exist.) Then we prove the proposition about solubility 
of equation p8|) . 



Lemma 3.1. The following asymptotic formula holds 

a(r) = |r|C/'(/i) + 0(l) as 



+ 00. 



(49) 



Proof. In order to prove this formula we consider the solution 7(y, t) of the boundary value 
problem and investigate its asymptotic behaviour as t — > +00. Note that 7(y, t) is a smooth 
function of both variables for large values of r, because the operator d^ /dY^ + oj'{U) has only a 
finite number of the Dirichlet eigenvalues. Let us write 



sinh Yt 



sinh hr 

and find the asymptotics of j{Y, r) for large |t|. Since this function satisfies the following problem: 

sinh Yt 



sinh Ht 



for re {0,h), 7(0) = 0, 7(/i) = 0, 



we have 



(7')^dr + r2 / ~2^Y 



uj'{U(Y))^^ dY 



J{U{Y)) 



sinh Yt 
sinh /it 



7dr . 



This is obtained by multiplying the equation by 7, integrating over (0, h) and using the boundary 
conditions after integration by parts. 

Let Cu) bounds cj' from above, then the last equality yields that 



h r 



(f) 



-7 



dY < 



7^^ + 



4C2 



sinh /it 



sinh^ yrdy, 



provided r is sufficiently large, and so the following estimate holds 

7(r,r) = 0(r-2) 



as r — >■ 00. 



(50) 



(51) 



Indeed, inequality (|50l) gives 



and so 



(7? + 



max f(Y,T) <2 \ [ 

yelo.h] \Jq 



dY < 



1/2 



1/2 



(7')'dr 



7^ dY 



< 



8C2 



for large values of |t|. The last inequality immediately yields ([51 
Furthermore, we have 



sinh Yt 
sinh /ir 



Jrl(Y-h) 



o 



(e-l-l") 
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which combined with ([?T|) gives that 

7(y,r) = el^K^^'') + 0{t-^) as |r| oo. 

Substituting this into the formula for a obtained in section 2.4 [its principal integral term is given 
by the second relation (|T7))]. we get 

cr(T) = — [ el^l(^-'^)r[/'(r)dr + 0(1) aslrl^cxj. 
h Jo 

The asymptotic formula p9)) is a direct consequence of this representation. 

Now we turn to the assertion about the behaviour of o'(r) when is close to a Dirichlet 
eigenvalue of the operator d^/dY^ + ljj'{U). Prior to formulating the result, we notice that every 
such eigenvalue is simple (provided there exists any). 

Lemma 3.2. Let the operator d^/dK^ +w'(?7) considered on (0, h) have as a non-zero Dirichlet 
eigenvalue. Then the following asymptotic formula holds 

^(^^ ^ _«[2<M! + 0(1) asT^n. (52) 
(t — j 

Here J*{Y) is the corresponding eigenfunction normalized in L^{0,h), and ^'^.{li) ^ 0. 

Proof. The constant 7^ {h) is not equal to zero because otherwise the function 7* has the zero 
Cauchy data aX Y = h, and so vanishes identically on (0, h) which is impossible. 

It is clear that is not an eigenvalue of d^ /dY^ + io'{U) provided r is sufficiently close to r*. 
Let 7(y, r) be a solution of problem (fT9|) for such a value of r. If we show that 

^(Y^^) = Z2:tihhl^+^(^Y,T) as r^r., (53) 

then ([5^ follows from this formula in view of the definition of a [see (ITS)) and (dH)]. Notice that 
the first term in the right-hand side is invariant under changing of the sign of 7* and the remainder 
v{Y^ t) is a smooth function of both variables. 

In order to prove the asymptotic formula (I53p we use the following representation 

l{y.r) = ^S^+v{Y^r) 

for r close to r*; here C is a non-zero constant. In order to find C we substitute the right-hand 
side into (IT51) and get that v must satisfy the following problem: 

~v" + [t^ -uj'{U)]v = ~C-i^, w(0)=0, v{h) = l. (54) 

Therefore, viY, r) exists only if C = — 7*(^). Indeed, multiplying ([54| by 7*, integrating the result 
over (0, h) and then putting r = r*, we obtain that 

C= -i,{v" -[tI -J{U)]v)dY = (w'S - 7"«) dr, 
Jo Jo 

because 7* is a normalized Dirichlet eigenfunction of d^/dY^ +10' [U) corresponding to r^. Now, 
we integrate by parts in the last integral and take into account the boundary conditions for v and 
7*. This immediately gives that C = —^i{h). Then formula ([53l) follows from the representation 
for 7 with the function v{Y, r) found from problem ([54]). It is clear that this function exists and is 
smooth. 
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Now we are in a position to prove the proposition about the solubihty of the dispersion equation 



Proposition 3.3. Let assumption (I) hold, and so there exists s.^, that solves equation (|14p for 
some 31. Let ([/, h) be the stream solution corresponding to s* . Then the following assertions are 
true for the dispersion equation with a defined by {U, h). 

(i) If the operator / dY^ +ui'{U) considered on (0, ft.) has no Dirichlet eigenvalues, then 



d3?, , /dh. 
ds / ds 



(55) 



Moreover, if 

^(..)/^(..)>0, (56) 
then equation (llSp has at least one positive root. 

(ii) Let the operator d^/dY^ + w'(?7) considered on (0, ft.) have no zero Dirichlet eigenvalue. If 
this operator has exactly k positive Dirichlet eigenvalues, then (jl8p has at least k positive roots. 
Moreover, if inequality (|56p holds, then (jlSp has at least fc + 1 positive roots. 

Proof, (i) Since the operator d^/dF^ +uj'{U) has no Dirichlet eigenvalues, cr is a smooth function. 
Therefore, it is sufficient to evaluate tT(0) in order to prove that guarantees the existence of a 
root. First, we consider the case when k = U'{h) > 0, and show that condition ([56]) is equivalent 
to the following inequality: 

cr(0) < 0. (57) 

Then assertion (i) is an immediate consequence of the last two inequalities and the asymptotic 
formula (|49p. according to which cr(r) +oo as t ^ +oo 
For proving ([55]) and ([57]) we note that 

cr(0) = K7'(ft,0) - +w(l) (58) 

by the definition of cr. On the one hand, 'y{Y, 0) satisfies the boundary value problem (fT9|) with 
r = and {U{Y; s*), ft(s*)). On the other hand, the Cauchy problem for a general stream function 
U{Y; s) is as follows: 

U" + uj{U)^0, C/(0) = 0, C/'(0) = s. 
Differentiating the first two relations with respect to s, we get 

[/" + a;'(;7)C/ = for C/e (0,ft), f/(0) = 0. 

(Here and below U denotes the derivative with respect to the parameter s.) Comparing these 
relations and the problem for j{Y, 0), we see that 

^(r,0) = 3^^^, (59) 
U{h{s^); s*) 

where the denominator does not vanish because jiY, 0) is well-defined by the boundary value 
problem with t = 0. Now we differentiate (15^1) and substitute the result into where k is 
changed to J7'(ft(s*)) in the first term. Then we get that 



^^P^ ^ U'{h{s.))U'{h{s^):s,) _ 1 ^ ^^^^ ^ [{U'(Ms))r] 



U{h{s.,);s,) K 2U{h{s);s) 



i+.(l). 
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Applying relation ([S]), we obtain 



a(0) 



1 



coil) 



U{h{s^.); s*) 

Here U'{h{s^,)) is changed back to k; moreover the equality /i(s*) = 1 and the formula 

U'ihis,)) h{s,) + Uihi.s,);s,) ^ 



(60) 



are used. The last equality arises when one differentiates the condition U{h{s); s) — 1 (it holds for 
every stream solution) with respect to s and puts s = s*. Finally, differentiating with respect 
to s (see formula dH]), where H the is defined), we get that 



1 



2 

3k 



d^ 
d7 



(..)/S<-) 



(61) 



In view of this formula, inequality ((57|) is a consequence of ((56)) when k — U'{h) > 0. This proves 
(|57p and completes the proof of the proposition when the last inequality is assumed to hold. 

It remains to consider the case when k = U'{h) < 0. According to the asymptotic formula (|49p . 
this condition yields that cr(T) — >■ — oo as r — > +oo. Hence, if wc show that inequality (1561) implies 
that 

f7(0) > 0, (62) 

then we immediately obtain the required assertion. Using formula (j6ip (it is independent of the 
sign of k) and the inequality k < 0, we see that ([62)) is again equivalent to ()56)) . Thus the proof of 
assertion (i) is complete. 

(ii) First, let exactly one Dirichlet eigenvalue > of (P/dY^ +Ld'{U) exist. Then formula 
(|53)) and the definition of a yield that (t(t) — > =Foo as t — > + provided ±k > 0. On the other 
hand, cr(r) — >■ ±oo as r — > +oo by virtue of the asymptotic formula ()49p . By the assumption there 
is no eigenvalue other then r^. Hence a{T) is smooth for t £ (r*,+oo) and it tends to opposite 
infinities at the ends of this interval. Then at least one root of equation ()18p exists on (r*, +oo). 

Since zero is not an eigenvalue, cr(0) is defined. Therefore, it is easy to modify the proof of (i) 
so that it will combine inequality ()56l) and the limit of ^{t) as t — >■ r.^ — instead of the limit as 
T — >■ +00. On this way, one obtains that one more root of (jlSp exists on the interval (0,r*). 



Now, let us assume that there are exactly two Dirichlet eigenvalues [t* ] and [t* ] 
Let ji^^ and 7!^^ respectively, denote the corresponding eigenfunctions. Then it is easy to see that 



air) 



{li'^Yih) 



0(1) 



as r 



JJ) 



i = l,2. 



which is similar to (15 



of k) as r goes to ri"^^ + and ri^' 
the interval 



Thus, cr(T) tends to opposite infinities (their signs depend on the sign 
0. Hence there exists at least one root of equation (|T8)) on 



J2) 



[ri^^T^'^^) . As in the case of a single eigenvalue, at least one more root belongs to 
-00) , and so the total number of roots is not less than two. The case of k eigenvalues should 



be treated in the same way. The proof is complete. 



3.2 The dispersion equation for unidirectional flows 

First, a stream solution ([/, h) describes a unidirectional flow, that is, U'{Y) > for all Y G [0, h], 
only when this solution depends on s > sq which is a root of equation ([T4| with 3?q^'(s) (see [50] . 
sections 5.2 and 5.3). Moreover, we have (see jSOj, section 5.1): 



(+) 



ds 



(s) > provided ± (s 



> 0. 



(63) 
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Second, the corresponding operator d^/dF^ + io'iU) considered on (0, h) has the empty Dirichlet 
spectrum. Hence the solution ^{Y] r) of problem ([T9| is defined for all r € R and is a smooth 
function of both variables. The same is true for a that stands in the dispersion equation because 

a(T) = a(0) + K v'{h; t), where v{Y; r) = 7(^5 r) - ^{Y; 0) - ^{Y; r) - ^ , (64) 

U{h) 

where the last equality is a consequence of formula (j59l) . Of course, v depends also on the parameter 
s > So and, for the sake of brevity, this dependence is indicated by the dot standing for the 
differentiation with respect to s. Note that this function does not vanish identically on (0, h). 

Lemma 3.4. Let s > sq he a root of equation (I14[) with Sig^' in the left-hand side. If the function 
V is defined by formula (|64p . then v(Y;t) < for all Y S [0,h] and v'{h;T) > 0; both inequalities 
hold for all t G H. 

Proof. The definitions of U and 7 and formula (|59l) yield that v satisfies the following problem: 

v" + uj'{U) V = T^-f on (0, h), v{0) = v{h) = 0. 

Note that 7 > on [0, h], which is a consequence of formula (3.9) in [5D]. Putting w = v/U' and 
taking into account the equality U'" = —U'uj'{U), we rewrite the previous boundary value problem 
as follows: ^ 

w" + 2^w' on {0,h), w{0)=w{h)=0. 

The right-hand side of the equation is non-negative because U'{Y) > for all Y G [0, h]. Therefore, 
the maximum principle is applicable (see [53], Chapter 1, Corollary to Theorem 4). Hence w < 
on (0, h) and the same is true for v. Combining this and the boundary conditions for we get 
that v'{h;T) > for all r e M. The proof is complete. 

Now we are in a position to prove the result concerning the dispersion equation (jlSp for unidi- 
rectional flows. 

Proposition 3.5. Let s > sq be a root of the equation 5l^\s) = r for some r > r^.. If the function 
a is defined by the stream solution {U, h) corresponding to this s, then the following assertions hold. 

(i) For any s € {so,Sc) equation (jlSp has a positive solution. 

(ii) Equation (jlSp has no positive solutions for s > s^- 

Proof. Assertion (i) (it is included for the sake of completeness) is an immediate consequence of 
Proposition 3.3, (i). Indeed, inequality ([55| follows from with the lower sign and 

^(s) < 0. (65) 

The last inequality holds for all s > sq according to formula (4.1) in [50] . 

In order to prove (ii), let us take an arbitrary s > Sc and assume that the corresponding 
equation (|T8]l has a positive solution tq. Writing formula (|48]) in the form 

ct(to) = -w(l) + TTT - Ky{h{s),To) = -uj{l) + - Kj'{h{s),0) ~ Kv'{h{s),To), 

n[s) n[s) 

we get that 

(T(To)<-c^(l) + -^-/t7'(/i(s),0) (66) 

because k — U'{h{s);s) > and Lemma 3.4 gives that v' {h{s);To) > 0. Using formula (|59p . we 
transform the last term in (1661) as follows: 



-.7'(^0) ^ -U'ik)Zil^ = . _s-.iumuiH) 

^ ' ' ' U{h) 2U{h) U{h) 
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Here the last equality is a consequence of equation ([5]) . Substituting this into (|66p and taking into 
account that U{h) ~ 1, we get that 

. , Kj S Ki S 

which (|60|) is also used. Subtracting equation (fT8| with r = tq from this inequality, we obtain that 
1 + (s/h) > 0, but this cannot be true. Indeed, it is shown in the proof of Proposition 3.3 that this 
inequality is equivalent to 



-(.)/^^W>». 







ds 

However the denominator of the last inequality is positive by ()63p . where the upper sign must 
be taken, whereas the numerator is negative by (|55|). Thus, the obtained contradiction proves 
assertion (ii). 

4 Proof of Main Theorem 

Our proof of Main Theorem is based on the following theorem which deals with bifurcation from 
a simple eigenvalue. 

Theorem 4.1. Let I be an open interval o/M such that £ /, and let X, y be Banach spaces. If 
a continuous map If : / x X M> ^ has the following properties: 

(i) the equality 9^(A, 0) = holds for all X G I, 

(ii) the operators 3^x, iF^ and 9^ac exist in a neighbourhood o/(0,0) and are continuous there, 

(iii) 3^(^ is a Fredholm operator, the null-space 0/3^^(0,0) is one-dimensional, 

(iv) if the null-space o/ 9^^(0,0) is generated by C,'^^\ then C}'^^ does not belong to the 
range of 3'^(0, 0). 

Then a sufficiently small £ > exists and a continuous curve 

{{\{t), at)) : |t| < e} C / X X, 
bifurcates from (0, 0). Moreover, for pairs belonging to this curve the following properties hold: 

^(t) = t(;(0) +o(i) whenO<\t\<e, 

and 

{(A, 0(^V:C^O and ?(A, C) - 0} = {(A(i), CW) : < |i| < e}, 
where V d I x X is a certain neighbourhood of (0, 0). 

This result is proved by Crandall and Rabinowitz (see Theorem 1.7 in |21j). 

4.1 Application of Theorem 4.1 to equations (1371) and (I38p 

In order to prove the first assertion of Main Theorem, we apply Theorem 4.1 to equations (|37l) and 
([55)1 . (It is sufficient to consider only the first of them.) 

Let us denote the left-hand side of equation (|37|) by F{X,C,). This is an operator mapping 
IR X H^'" I— H]y'". Indeed, the derivatives of (f> in the expression for F depend on C and A through 
the weak formulation ([35]). According to Remark 2.3, assumption (III) implies that <p G Hq'^^(S') 
provided C is sufficiently small. Here Hq'^^(5) denotes the subspace of H^'"(S') whose elements 
vanish on dS, that is, Hq'^^(S') consists of C'^'^-functions on S that are Ao-periodic and even in 
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X, and vanish for z — and z = h. Thus, conditions (i)-(iv) of Theorem 4.1 must be verified for 
the equation F{\, () ~ 0, in which case X = 11^'" and y = njj". 

Condition (i) of Theorem 4.1 foUows from assumption (III) because 

F(A, 0) = h) for A e /, x G (-Ao/2, Ao/2), (67) 

and we have that (j)z{x,h) — when — 0. In order to prove this fact, we substitute these two 
equahties into the integral identity ([34|) . thus getting 

pAo/2 

/ / + (1 + >^f(t)zVz + (1 + \f Mu) ~uj{u + (j))]v) Axdy = 0, 

J -Ao/2 Jo 

where v belongs to the same class as in ([M)) . Since lu E C^'°'(R), the integrand can be written as 
follows: 

(j)^ + (l),Vz - Lo'{u)(j)v + (2 A + A^) [(j),vz - Lo'{u)(j)v] + (1 + \fv 0{(f>^). 

Let £ denote the operator of the Ao-periodic boundary value problem (l43l) (note that L corresponds 
to the first three terms in the last expression). Then assumption (III) implies that 

L : IIp'^^ (5) I— > IIp'^^ (5) is an isomorphism 

(see, for example, [1]). Therefore, we have that 

ll<^lln^'° IS) — C'||'/>||n2,a r-SM where C = const independent of (f>. 

Hence </> is identically equal to zero provided is small. From this fact and (j67| one gets that 
condition (i) is fulfilled for F. 

Condition (ii) concerns the existence and continuity of the operators F\, F(^ and Fx^ in H^'". 
According to Remark 2.3, assumption (III) implies that the mapping (A, () cj) [it is defined by 
the weak formulation psp ] is continuous from a neig hbourhood of (0, 0) in IR X H^;" to U^'l^iS). 
Moreover, this mapping has continuous derivatives of the first and second orders. 

Condition (iii). The fact that F(^ is a Fredholm operator follows from (|37|) . Furthermore, 
Proposition 2.2 yields that 

[Fc{0,0),]{x)^(x^^k){x)-(j-~^y{x), xe (-Ao/2, Ao/2), (68) 

where (%^^\) (x) is defined as follows. Let $(^)(x, z; <;) be Ao-periodic in x and satisfy the integral 
identity (gHl) with C = ^, then 

{x)=^i^\x,h;,), 

and so X*^^^ is a Fredholm operator from H^'" to 11^^'". 

Finally, conditions (a) and (b) in assumption (II) guarantee that the null-space of F(^(0,0) is 
one-dimensional. Indeed, it is straightforward to check that this space is generated by 

C(o)(,)^eos^. 

Ao 

Condition (iv). In order to obtain an expression for the operator Fx^(0,0) let us consider ^p^'^\ 
the second-order term in the expansion of 0. It follows from (1551) that ip'-^^ satisfies the following 
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integral identity: 

fAo/2 fh 



(2) , 



-Ao/2 JQ 
Ao/2 



Ax dz 



cj(m) 



Ao/2 Jo 

■Ao/2 fh 



2AC+ (2A- J )C 



dx dz 



w > dx di 



Combining the last equality, the second expression (pi]) and formula we get that 
[F,c(0, 0) ^] (x) = (X^'k) (x) - ( J - ^) ^(^), X e (-Ao/2, Ao/2). 



(69) 



Here (X^^\) (x) is defined in the same way as the operator X^^K Namely, let $^^''(x, z; <^) be a 
Ao-periodic in x solution of the problem that has the following weak formulation: 



Ao/2 fh 



Ao/2 Jo 

■Ao/2 /.h 



$^2) + $(2) _ $(2) „ 



da: dz 



Ao/2 Jo 



+ ^i^\x, z; <r) - ^i^^x, z; ^ + uj' {u)<^>'^^\x, z; w 



da; dz . 



Here v belongs to the same space as in (|34p and ^^^-^ (•;(;) is defined after formula ((68| . Then we 
have that 

(X^^\){x) = '^f\x,h-c,). (70) 

Now we are in a position to check condition (iv). Since is a Fredholm operator, it is sufficient 
to show that 

W(o,o)c("\C^°))^o, (71) 



where (•, •) denotes the corresponding bilinear form. In fact, (j7ip is a consequence of condition (a) 
in assumption (H) as the following assertion demonstrates. 

Proposition 4.2. In formula (j7ip . the left-hand side is equal to 7ro''(To). 

Proof. It follows from ^ and dTO]) that 



because 



Ao/2 



cos 



o 27ra: , Ao 
da; = — , 



Ao/2 



Ao 



Seeking {x,h;C'-°^), it is natural to use the same ansatz for both and namely: 



fx,z;C(°)) =p^'\z) cos^, 
V / Aq 



Then we get that 



(^^Ac(o,o))d°)],C(°))-^ 



fc = 1,2. 



K, 1 
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where p^^-* satisfies the boundary value problem: 

Furthermore, the problem for p^^^ is as follows: 
which, in particular, implies that 

p«(/i) = -ao(ro) 

[see the derivation of formula (|47|) ]. 

Putting p(^) = 2^*^^' + g and taking into account the last equality, we see that 



Fac(o,o)c(°),C(°)) = ^ 



K 1 

h K 



(72) 
(73) 

(74) 



where [cf. formula (|T7|) ] 
because g satisfies the following problem: 



-q,, + [t^ - ij'iu)] q = + [t2 + uj'iu)] p^^^ - - t^zu,, q{Q) = q{h) = 0. 

Splitting the last integral into the sum of two terms and integrating by parts twice in the first one, 
we obtain 

q,(h) = - f [-i,,{z,To)+J{u{z))-f{z,To)] p«(^)d^ 

Jo 

2 



To / 7(z,To) 







P^^\z) - -ZUz{z) 



Az. 



In view of equation for 7 [see (flQ)) ]. this reduces to 



q,{h)^~p(^\h)^2Tl j{z,To) p^'Hz)-h-'zu,iz) 
Jo L 



dz, 



which combined with formulae ((73)) . ((74|) and the definition of tq leads to 

(f,c(0,0)C(0), C(°)) = -AoT^ r7(z,ro) [pW(z) - /^-lz^i,(z) 



dz. 



(75) 



Let us show that the last integral is equal to — o-q(to)/ (2 tq), which completes the proof in view of 
the relation between tq and Aq, and the equality — cro(r) = (t'{t) [see formula (|48l) ]. Differentiating 
formula ((47| . we get 

'^o('^o) = ^ W {27(z,to)tozm2(z) +7'(2;,ro)[roZu^(z) + 2a;(?i(z))]}d2:, 
Jo 

where 7' stands for d'y/dr. Let us also differentiate relations ([T9|) with respect to r, thus obtaining 
-7L + k'-^'H]7' = -2r7, 7'(0) = 0, y{h)^0. (76) 
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This problem is similar to (|72l) . and it allows us to express 7' in terms of 7 (of. section 2.4): 



-2r 



G{z, z'; t) 7(2;', r) dz', 



where G(z, z'; r) is Green's function of problem (1751) . Using this representation for calculating the 
second term in braces in the formula for (To(ro), we get that this term is equal to 

—r^ [tqZU:,{z) + 2uj{u{z))]dz G{z,z';T)j{z',To)dz'. 



After changing the order of integration, we notice that the inner integral gives hp^^\z'). Thus, 
the last double integral reduces to 



2 TO / p(i)(z')7(^',ro)dz'. 
Jo 

Substituting this into the formula for (t'(to), we arrive at 

c^o('^o) = 2ro / l{z,To) p^^\z) - h^^zu:,{z) 



dz. 

Comparing this and (|75l) . we obtain the proposition's assertion. 

Since conditions (i)-(iv) of Theorem 4.1 are fulfilled for equations ((37|) and ((38)) . assertion (i) 
of Main Theorem is proved. 



4.2 Proof of Main Theorem: assertion (ii) 

It is noticed in Remark 2.3 that there exists a solution of the nonlinear problem whose weak 
formulation is given by the integral identity (j35p . Moreover, the map {\{t)X{x,t)) 1— > (t){x,z,t) 
defined by (|35p is continuous from a neighbourhood of (0,0) in IR x H^'" to the space ng'^^(S'), 
and the following estimate mentioned in Remark 2.3 holds: 

Aq ^ ' Aq 

Now we change variables using formula (1^^ . where A = Ao[l + A(t)], and (PT|) and p2p . This gives 
the function 

*<^- « ^ " + * {^■'^mry') " ""<*'^ 

In view of considerations presented in sections 2 and 4.1, ^I^ solves problem ([l])-(l3]), whereas the 
pair (^',^) satisfies condition (U). 

Comparing formulae (|77l) and (|22p . we see that the second term in the right-hand side of ([Tf)) 
is equal to in ((22|) . Combining this fact, the last inequality and assertion (i) of Main Theorem, 
we get that the second term in the right-hand side of is 0{t). This completes the proof of 
Main Theorem. 



5 Examples 

In this section, we consider two examples of vorticity distributions (the corresponding flows of 
constant depth were studied in our paper |50|). for which the dispersion equation is investigated 
along with condition (l56l) . This allows us to draw conclusions about the existence of Stokes waves 
perturbing a flow of constant depth with vorticity. In the first example, the vorticity is equal 
to a positive constant, while the second example deals with the linear vorticity having a positive 
coefficient and vanishing at zero. 
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Figure 1: For the constant vorticity uj — b > 0, the functions "R^^^ (dashed hnes) and D?^^' (bold 
hnes) are plotted for b = 1/2 (left) and b — 2 (right). 



5.1 Constant positive vorticity 

Let UJ = bhe a positive constant, then it immediately follows (see [50] section 6.1) that so = V^b, 
and 

s 



4+)(s) 



2b 



+ Vs 



2b 



(78) 



are the only non- vanishing depth functions defined for s > sq and such that h'^\s) < h\^\s) 



with the equality holding only for s = sq when both values are equal to ho = \/2/6. Hence the 
left-hand sides in the Bernoulli's equation are as follows (see Fig. 1, where these functions are 
plotted for 6 = 1/2 and b — 2; reproduced from [SD], p. 391): 



2 + 2 6 



These functions have the following properties: 31q'''-'(s) < 3?^"^''(s) for s > sg, while 



T - 2b 
b 



(79) 



(80) 



and 



2 "T 2 



2 

y + ^(sTs) + 0(l) ass- 



-l-oo. 



(81) 



The distance between 3lo"*'^(s) and R['^\s) increases with s due to the second term in the right-hand 
side. [This term is lost in formula (6.2), [50], analogous to ((8T]) .] Furthermore, we have that 



2 ~r 2 

ds 



(s) — =Foo as s — )■ So + and 



(+) 



,i — (s) > for s > So. 
ds 



(82) 



The first of these relations shows that the graphs of JIq^'' and have the common vertical 

tangent at the point (so,ro), while the latter inequality implies that the function D?^^"* increases 
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strictly monotonically, and so its graph of lies strictly above the level tq. Finally, the single 
minimum Tc < ro of J^q^'* is attained at Sc > sq satisfying the following equation: 

is.br-2b'=(-^X . (83) 



1 + Scb 

It is clear that the first component of the stream solution is 

U{Y;s) = -^Y^ + sY . (84) 

It solves the corresponding Cauchy problem [see ([6]) for the Cauchy data]. Since there are two 
options for the second component [they are given by formulae (1781) ]. two steady flows of constant 
depth exist for every r > r^. It is shown in ISO], p. 392, that if the depth is equal to /iq^-'(s) and 
s € (sq, Sc) is obtained from equation (|14p with 3?o^'' [it will be clear from what follows why s > Sc 
is not considered], then the flow is unidirectional because U'{Y; s) > for all Y € [0, h'^\s)]. On 
the contrary, if the depth is equal to h'{^\s) [in this case, s > sq satisfies equation ((Ti)) with 3^1^''], 
then there is a near-surface counter-current. Indeed, it is easy to calculate that 

U'{h'i^\sy, s) = - Vs^ - 26 < and [/'(F; s) < for F > s/& 

It is clear that the latter value is less than h[^\s). 

In order to write down the dispersion equation we find the following solution of problem ()19p : 

sinhrF , , , ^ 

Here and below h — h{s) is either hf^\s) or h[^\s) [see ((78|) ]. and s is a root of the corresponding 
equation (IHl) : the common form of both these equations is as follows: 

s'^ + 2h{s) = 3r + 2b. 

We substitute k — —{bh — s) and 7 into (fT8|) and obtain 

T{bh - s) cothr/i - 6 - [bh - s)-^ = 0. (85) 

This dispersion equation coincides with ((T7)) when 6 = 0. Indeed, we have that s = h = /iq^'' in 
view of the first formula ((78)) [the second formula (|78l) is meaningless for 6 = 0]. 

It is easy to see that the derivative of the left-hand side in (1551) has a definite sign that depends 
on the sign of bh — s, and so if this equation has a root, then it is simple. Instead of demonstrating 
directly that ([85]) is soluble for both hf^\s) and h'l^\s) [see ([78)) ]. we apply Proposition 3.3 for 
this purpose. 

First, let r > tq, then the second relation (|82p implies that the function 3?^^-* increases mono- 
tonically. Thus, the inequality 

d/ii+) , . 1 



'—(s) = ± 1 + ^ " >0 
us ^ ^ 6 V V^^^2bJ 

guarantees that condition (|5B| is fulfilled, and so Proposition 3.3 yields that ([55)1 with h — h'{^\s) 
has a positive root. In this case, the left-hand side in (|85p is a monotonically increasing function 
of T because bh — s = \/s^~~2b > 0. Hence the dispersion equation (|85p with h = /i^^'' (s) has only 
one positive root. 
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Let us turn to the case when h = hQ'^\s) and s G (sq, Sc). Then we have that r G (''c, ?'o) and 

d4+) , . 1 



(s) = - 1 -^=^= < for all s > sq. 



ds ^ ' b \ Vs^ - 2b y 
Moreover, in view of ([55| the following inequality holds 



5 + T ( 1 



d3?f+) 



ds ^ ' 3 



< for S e (sq, Sc 



Hence Proposition 3.3 yields that the dispersion equation ([55]) with h = hQ^\s) has a positive 
root. Again, there is only one root because the left-hand side in ([55)1 is a monotonically decreasing 
function of r, since bh — s = —y/s^ — 2b < in the present case. 

On the basis of the above considerations we formulate the following. 

Proposition 5.1. Let lo ^ b, where b is a positive constant, then the following assertions hold. 

(i) For every r > tq [see (j80p for ro] the equation 3i\^\s) = r [see (|79p for 3?^^'] has the unique 
root sl^'' > Sq. In the flow of the constant depth h — li^^^ ('^*^'') ^^^^ second formula (|78p], the 
velocity field is described by the stream function The dispersion equation with s = s^^^ 
and h = h^^"' (**^^) ^'^^ onZy one positive simple root. 

(ii) For every r S (fc'^o) [relations (j83p and (|79p </iat give Sc and Jt'^K respectively, show 
how to find rd the equation 5l^\s) = r has the unique root sl^"* G (so,Sc). /n </ie /Zow of the 
constant depth h = /iq^' (^**'^) ^^^^ /irsf formula (|78p ] . t/ie velocity field is described by the 

stream function (|84|) . T/ie dispersion equation ([85]) mi/i s = sl^"* anrf /i = /iq^'' (^s*°'') ^'^^ only one 
positive simple root. 

Combining this proposition and Main Theorem, we arrive at the following conclusion. 

Corollary 5.2. Let u ^ b, where b is a positive constant. Then for every r > r^, r ^ r^, there 
exists a family of Stokes waves with wavelengths belonging to a neighbourhood of Aq such that 
27r/Ao is the only positive root of equation ([55]) depending on r as described in Proposition 5.1. 
These small- amplitude waves perturb each of the flows whose constant depths are given by formulae 
(j78p . where s depends on r through the corresponding equation (|14p . 

Remark 5.3. When w = 6 is a positive constant, the family of Stokes waves, existing for every 
r > ro by Corollary 5.2, perturbs a flow of constant depth that has a near-surface counter-current. 
If r G (rc,ro), then the corresponding family of Stokes waves perturbs a unidirectional flow of 
constant depth. In the case when uj = b < 0, our approach gives the results similar to those 
formulated in Proposition 5.1 and Corollary 5.2. The corresponding flows of constant depth (they 
have near-bottom counter-currents), are described in f50J, section 6.2. 

Remark 5.4. Further details about small-amplitude Stokes waves with constant vorticity can 
be found in ^68^ and [20j. In the paper (68j, Wahlen paid much attention to the behaviour of 
streamlines for such a flow. In particular, he investigated in detail the streamlines that form a 
'critical layer' which separates two layers with the opposite directions of flow. Moreover, within 
this layer, all streamlines are closed, that is, it consists of the so-called cat's-eye vortices (see 
figures 1 and 3 in [68| . where w = 6 < 0). Besides, a full description of particle paths is given in 
[68] . The approach developed for the same problem by Constantin and Varvaruca [20^ is applicable 
to a wider set of free surface profiles, in particular, it includes overhanging ones. 
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Figure 2: For the vorticity w = &t, & > 0, the interlacing pairs of functions 3^2fc (dashed hnes) 
and 3?2^^]^ (bold lines), and 3^2*-^ (dashed lines) and 3^2fe+i (bold lines), A; = 0, 1, ... , are plotted 
for 6=1/2 (tips are closer to the vertical axis) and 6 = 2. 



5.2 Linear positive vorticity 

Let uj = 6t, where 6 is a positive constant. Then VI{t) — 6r^/2, and so sq = ^/h. According to 
formulae (4.4), (4.5) and (4.8) in [50 , the following two sequences 

/i(+'(s) = ^-|^arcsin^+j^ and hY\s) = hf\s) + J=0,1,..., (86) 

give the depths of flows with horizontal free surfaces. Therefore, the left-hand-side terms in equa- 
tion p4)) are as follows (see Fig. 2, where the corresponding graphs are plotted for 6 = 1/2 and 
6 = 2; reproduced from [50)): 



1 . 2(-l)^' 

s —o^ ;= — arcsm 

V6 



yfb 2-K ( 11 
V + 7^(,^'+2^2 



, i = o,i. 



(87) 



In particular, we have that 

= i - 6 + ^ arcsin ^ j and 3?i+^(s) 



2 . \/6 
^ —= arcsm ■ 



, 27r _ ... 
6 H — — arcsm 



(In each of the two families shown in Fig. 2, the graphs corresponding to these functions are plotted 
as the lowest dashed and solid lines, respectively.) Hence 



3%/6 



whereas Tc — the only minimum of Jl^\s) — is attained at 



\2 
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Furthermore, both relations ([5^ remain true in the present case together with the conclusions 
drawn from them about the pattern formed by the graphs of 3?o~'''* and Jl^i'^ . Moreover, this pattern 
is repeated infinitely many times; indeed, formulae (j87p imply that it is shifted by 27rj/(3-\/5), 
j = 1, 2, . . . , along the positive r-axis (see Fig. 2). 

The first component of a stream solution has the form 

U{Y;s) ^ ±^ sin VbY . (89) 

Indeed, this function solves the Cauchy problem that consists of the first and second relations ([5]). 
The sign in (1891) coincide with that in the superscript of the corresponding /ij^^(s), which is the 
second component of the same stream solution. Infinitely many options of it are given by formulae 
([86|): however, only a finite number of them defines steady fiows of constant depth for a particular 
value of r greater than re- The number depends on how many roots has the corresponding set of 
equations (fT4|) for a chosen r. Note that the number of roots increases with r. 



Prior to investigating which fiows support Stokes waves we turn to the sequence of dispersion 
equations which exists in the present case. Indeed, the left-hand side of (fT8| involves 7 that solves 
the following problem: 

-7" + (r2-6)7 = 0, 7(0) = 0, -/(h) = 1, (90) 

where h — h'j^\s), j — 1, 2, ... . Each of these values depends on a root s of equation (fT4|l with 

31^*'' in the left-hand side. Applying the remark that follows problem ([T^. we conclude that ([50]) 
has a solution for all r such that 

rV6-(^)', fc = l,2,.... (91) 

This condition is obviously fulfilled when > b, in which case 7(Y,t) is equal to 

• Y/h when = b; 



sinh \/t'^ — bY o 

• j^=^= — when T > b, 

sinh \/t'^ — bh 



and so 7 (and cr as well) is a continuous function of t when > b 
we get that the sequence of dispersion equations for >b has the form: 

sVt2 - b cof? [ybhf\s)^ coth [^/t"^ -bhf\s)^ - 1 ± 6s cos [Vbhf\s) 



Using formula with h = h'f^\s) and the second of the buUeted expressions in formula psp . 



= 0, 



J =0,1,..., (92) 



and the left-hand side must be understood as the corresponding limit for = b. 
If < 6 and condition (IQTI) is fulfilled, then 



^ sin Vb-r'^Y 
sm y — T'^ a 

Substituting (j89p and this expression into formula (1181) , wc get the sequence of dispersion equations 
for this case: 



Vb - r2 cos^ (^Vbh^f'\s)'^ cot (^Vb~^ h^f'\s)j -l±bs cos (Vb/iJ^^^ (s)) =0, 

i = 0,l,... (93) 
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The hmit form of ([M)) as — ^ 6 — coincides with that of (IMI) as ^ 6 + 0. It is clear that the 
derivatives of the left-hand sides of (j92|) and (|93l) have definite signs on the intervals of continuity, 
and so if these equations have roots, then they are simple. In order to investigate equations (I92p 
and dnS]) we apply Proposition 3.3. 

Taking into account formula (|9ip . we begin with flows such that ^/bh'j^^s) < vr. Formulae 

((86)) show that only hQ^\s) and h\^\s) satisfy this inequality. Moreover, Proposition 3.3 (i) 
guarantees that for each of these depths equations (|92p and have at least one positive root 
provided condition (|56p is fulfilled for either of the following pairs: 

= \ h = Hq ^5 3^ — "J^^ \ h = /i^ ^ • 

(In fact, the form of the left-hand side in and allows the existence of only one root.) The 
results obtained in [50], section 6.3, yield that inequality ([56]) always holds for the second pair of 
functions, whereas for the first pair it holds only when r £ (rc,ro). [We recall that vq — n/Sy/b, 
whereas Vc — 31q'^\sc) and Sc is given by formula ([55]) .] 

In the first case, the depth is ft.Q^-'(s), where s £ {so,Sc) is obtained from equation (|14l) with 
IKq"'''* and r G (rc,ro). Any such flow is unidirectional because of the upper sign in formula ()89p 
which gives that 

U'{Y; s) = s cos VbY > for all Y e [O, hl^\s)] . 

On the other hand. Proposition 3.5 is applicable when 31 — Jl^^^ and h — hl^\ Therefore, if a is 
defined by the stream solution corresponding to the root s > of 3iQ^\s) — r with an arbitrary 
r > Tc, then the dispersion equations (|92t and (|93)) have no positive solutions. 

In the second case, the depth is h^i^\s), where s > sq is obtained from equation with 
and an arbitrary r > tq. All such flows have a near-surface counter-current because 

U'{Y; s) = scosv^F changes sign only once on (O, h[^\s)). 

Combining the above considerations and Main Theorem, we formulate the following. 
Proposition 5.5. Let lo = br, where b is a positive constant, then the following assertions hold. 

(i) For every r £ (rc,ro) there exists a family of Stokes waves with wavelengths belonging to 
a neighbourhood of Aq such that 27r/Ao is the only positive simple root of (j92p and (j93p in which 
s £ (so, Sc) solves the equation 31q^\s) = r. These small- amplitude waves perturb the unidirectional 
flow having the constant depth h — hf^\s) [see the first formula (j86p ]. 

Equations (j92p and (j93p have no roots when s > Sc is a solution of the equation 31q"''^(s) = r. 

(ii) For every r > tq there exists a family of Stokes waves with wavelengths belonging to a 
neighbourhood of Ag such that 27r/Ao is the only positive simple root of (j92p and (|93p in which 
s > sq solves the equation 3l[^\s) = r. These small- amplitude waves perturb the flow with a 
near-surface counter- current; the constant depth of this flow is equal to h ^ /i[~*^'(s) [see the first 
formula (|86p ] . 

If Vbh'^'^\s) £ (7r,27r), then formulae (|5S|) show that only \s) and h[ \s) satisfy this 
condition provided s is a root of 3^q~^(s) = r and 3i^~''(s) — r, respectively. It is clear that the 
second of these equations has only one solution for any r > r^ -\- and the range of solutions is 

(sq, -I-oo). Furthermore, the equation !Kq \s) — r has only one solution belonging to (sc, -foo) for 
any r > Tc + Moreover, for every r £ (re + ■j^^j'^o + ^be latter equation has another 

solution in (so, Sc). Thus, the pair of equations 3?g \s) — r and '3l[ \s) — r has a pair of solutions 
for all r > Tc + -^jf except for r = tq -I- , when only the first of these equations has a solution. 
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All flows, whose constant depth is equal either to /iq (s) or to h\ (s), have near-bottom counter- 
currents. In the first case, there is no other counter-current, but the near-surface counter-current 
is also present in the second case. 

Now we substitute hg (s) and h\ (s) corresponding to the values of s obtained in the previous 
paragraph into equations (|92l) and (|93l) and see that Proposition 3.3 (ii) is applicable. Indeed, 
formula (|9T|) yields that = b— (tt/Zi)^ with h = h\. \s), fc = 0, 1, is the single Dirichlet eigenvalue 
of the operator d^/dF^ -I- (j-i'iU) considered on (0,/i). It follows from the cited proposition that 
equations ((92|) and (|93|) with h[. \s), fc = 0, 1, have at least one positive root t'^^ e (t*, +oo). In 
fact, the form of the left-hand sides in (|92p and (|93| yields the existence of only one root there. 

On the other hand, if condition ((56|) is fulfilled for the following pairs 

3^ = \ h — Hq \ and 3i. = \ h = \ 

then Proposition 3.3 (i) is applicable. From the results obtained in [50], section 6.3, it follows 
that (|56l) is always fulfilled for the second pair of functions, whereas for the first pair it holds only 
when s € (so,Sc), that is, when r G (tc + 27r/3-\/6, ro -I- 27r/3-\/6) in equation (fT4|) with 31q •* in 
the left-hand side. Then Proposition 3.3 (i) guarantees that equations (1^ and ([Mt with h\r\s), 
fc = 0, 1, have at least one positive root t^'^ e (0,t*). (Again, the form of the left-hand sides in 
and yields the existence of only one root there.) 
Combining the above considerations and Main Theorem, we formulate the following. 

Proposition 5.6. Let lo — hr, where b is a positive constant. Let s > sq be an arbitrary solution 
of \s) = r, k = 0, 1, with r > rc + 2'K/sVb (there always exist two such solutions except when 
r = ro + 27r/3\/5, in which case only one solution exists). Then the following assertions hold. 

(i) Equations ()92|) and (|93p with h — h\. \s) have a positive simple root t'^^ > \b — (tt/Zi)^] 
There exists a family of Stokes waves with wavelengths belonging to a neighbourhood of 27r /t|'^'' , 
and these small- amplitude waves perturb the flow of the constant depth h\. \s). Moreover, a near- 
bottom counter- current is present in the flow when k — 0, 1, whereas for k = 1 a near-surface 
counter- current also exists. 

(ii) Lf s > Sq is a solution of 3?^ \s) = r with r > Tq + 27r/3\/6 , or if s G (sq, Sc) is a solution 
of 'Jl'f\s) = r with r G (rc + 27r/3\/fe,rc + 27r/3\/6), then equations ([92]) and with h = h^ir\s), 
fc = 0, 1, have a positive simple root < [b — (vr/Zi)^] ^^"^ . There exists a family of Stokes waves 
with wavelengths belonging to a neighbourhood of2Tr/T_ , and these small- amplitude waves perturb 
the flow of the constant depth h\, ^ (s) . Again, a near-bottom counter- current is present in the flow 
when fc = 0, 1, whereas for k = 1 a near-surface counter- current also exists. 

If Vbh^j^\s) e (27r, 37r), then formulae ([55)1 show that only /i2^''(s) and h^^^\s) satisfy this 

condition provided s is a root oi3i^2^\s) ~ r and lK3^''(s) = r, respectively. The same analysis shows 
that these equations have the same roots as in the previous case provided 2n/3\/b is added to r that 
stands in the right-hand side of 3?o \s) — r and '3l[ ' (s) = r, respectively. But now, formula (|9l]) 
yields the existence of two Dirichlet eigenvalues, namely, t^^ = b—{ii/h) and r^2 = b—{2TT/h) with 
h — h^i^\s), k — 2,3, for the operator d^/dF^ -I- u;'(f7) considered on (0, h). Hence Proposition 5.6 
can be easily reformulated for this case. In fact, the form of assertion (ii) remains the same, but 
there is no near-bottom counter-current and the number of counter-currents is different. Besides, 
assertion (i) should be modified in view that more roots of equations (|92l) and (^5)) exist. The 
required amendments immediately follow from Proposition 3.3 (ii), and are left to the reader. 

In conclusion, we summarize what was said above about Stokes waves on fiows with positive 
linear vorticity. First, if r > rc, then every root of each equation 3?j*''(s) — r [see formula (|87|) for 
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"Jij (s)], defines a stream solution from wliicli small-amplitude Stokes waves bifurcate, unless this 

root is equal to so or is greater than Sc and satisfies the equation 3?q''^''(s) — r. Second, the number 
of stream solutions from which small-amplitude Stokes waves bifurcate increases as r grows. On 
the one hand, for larger values of r more equations 3i'j^\s) = r give roots that can be used in 
the dispersion equations (|92l) and (|93|) . On the other hand, the latter equations have more roots 
defining bifurcation wavelengths because larger roots of iK^^-'(s) — r exist for greater values of r. 
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